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Abstract 

The CUORE Crystal Validation Runs (CCVRs) have been carried out since the end of 2008 at the Gran Sasso National Laboratories, 
in order to test the performances and the radiopurity of the Te02 crystals produced at SICCAS (Shanghai Institute of Ceramics, 
Chinese Academy of Sciences) for the CUORE experiment. In this work the results of the first 5 validation runs are presented. 
Results have been obtained for bulk contaminations and surface contaminations from several nuclides. An extrapolation to the 
CUORE background has been performed. 
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1. Introduction 

The production of the CUORE [ 1 ] crystals was appointed to 
SICCAS (Shanghai Institute of Ceramics, Chinese Academy of 
Sciences) and began in March 2008, with the synthesis of the 
TeC<2 powder. Given the goal of CUORE in terms of back- 
ground (<0.01 counts/keV/kg/y at the Q-value J2l), the radiop- 
urity of the TeC<2 crystals is a crucial issue. 

Radioactive contaminations may come from long-lived, nat- 
urally occurring isotopes, such as 238 U and 232 Th and their 
daughters and from cosmogenic activation of the detector mate- 
rials after their production. To minimize the influence of long- 
lived nuclei, great care is devoted to the selection of all materi- 
als and ancillaries used for the preparation of the detector. Sea 
level transport and underground storage of prepared crystals are 
necessary in order to minimize their cosmogenic activation. 

A dedicated protocol [3J is defined for the radio-purity re- 
lated quality control of the crystal production process, start- 
ing from metallic tellurium synthesis to the final processing 
of ready-to-use TeC<2 crystals. Radio-purity certification pro- 
cedures, involving ICP-MS (Inductively Coupled Plasma Mass 
Spectrometry) measurements, y spectroscopy with HPGe de- 
tectors and a spectroscopy with Surface Barriers Detectors 
(SBD), are applied in each production phase to test the above 
mentioned materials. 

At the same time, cryogenic tests are designed to test the 
ready-to-use TeC<2 crystals upon their arrival at LNGS. The 
crystal validation is performed through experimental runs, each 
called CCVR (CUORE Crystal Validation Run), in which 4 
crystals randomly chosen from a batch coming from SICCAS, 
are mounted in an setup similar to a CUORE single module and 
operated at cryogenic temperatures for a time period of several 
weeks in order to test the bolometric performance and the com- 
pliance of the crystals to the contract limits in terms of radio- 
purity (see Table[TJ. 



Isotope 


Allowed Contamination 




<3-10- 1J g/g 


m Th 


<3-10- 1J g/g 


m Pb 


< 1 10 b Bq/kg 




< 0.1 Bq/kg 



Table 1: Contamination limits for the ready-to-use Te02 crystals (J. 



In this work the results of the first 5 validation runs, carried 
out from the end of 2008 to the middle of 2010, are presented. 
A summary of all CCVRs detector operation is reported in Ta- 
ble|2] The 4 crystals tested in CCVR1 (two of them were again 
tested in CCVR2) were sent to LNGS by plane. This was neces- 
sary in order to ensure a fast response on the radioactivity level 
of the crystals. All the other crystals were transported by ship. 



Run 


Duration 


Livetime [d] 


CCVR1 


Dec 20, 2008 - Mar 9, 2009 


59.9 


CCVR2 


Jun 6, 2009 - Jun 30, 2009 


19.4 


CCVR3 


Nov 11,2009- Jan 4, 2010 


43.05 


CCVR4 


Mar 31, 2010 -May 17, 2010 


25.8 


CCVR5 


Aug 11, 2010 -Oct 1,2010 


30.3 



1 Presently at: Joint Research Center, Institute for Reference Materials and i 
Measurements, 2440 Geel - Belgium ( 
Presently at: Nikhef, 1098 XG Amsterdam - The Netherlands 
3 Presently at: Institut d'Astrophysique Spatial, 91045 Orsay - France 
4 Presently at: CEA / Saclay, 91191 Gif-sur-Yvette - France ' 
5 Deceased i 



Table 2: Summary of CCVR data taking period and live time. 

42 Radioactivity study was performed on the high energy re- 

43 gion of the spectrum (above 4000 keV), where the contribution 

44 of the a lines from uranium and thorium decay chains is ex- 

45 pected. Given the short range of a particles, their signature 

46 is a clear indication of a radioactive contamination within the 

47 crystals or on their surface. Results (mainly upper limits) are 

48 obtained for both bulk contaminations (Sec. [3]) and surface con- 

49 taminations (Sec. [7). For 210 Pb, for which the study of the a 
so lines was not possible, the lower energy portion of the spectrum 

51 has been used (Sec. [6). 

52 An extrapolation to the CUORE background is performed in 

53 Sec. [8] 

54 2. Experimental setup 

55 Each CCVR setup consists of an array of four crystals, ar- 
se ranged in a single floor which in a first approximation repre- 

57 sents the single CUORE module. The four 5x5x5 cm 3 crys- 

58 tals are enclosed in a pair of Oxygen Free High Conductivity 

59 (OFHC) copper frames which serve both as mechanical support 
eo and thermal bath. The two frames are connected to each other 

61 by four small columns, also made of copper. Both frames and 

62 copper are wrapped with several layers of polyethylene. The 

63 crystals are connected to the copper frames by small Teflon 

64 supports that represent the weak thermal conductance versus 

65 the heat sink. 

ee Each crystal is provided with a Neutron Transmutation 
67 Doped Ge thermistor (NTD), which converts the phonon signal 
ea into a detectable voltage pulse |4|]. Some of the CCVR crystals 

69 are provided with two NTD thermistors. In these cases both 

70 channels are analyzed and the best performing channel (from 

71 the point of view of the energy resolution) is used for the final 

72 results. 

73 The CCVR experimental setup is hosted in a dilution refrig- 
erator placed in the Hall C of National Laboratory of INFN at 
Gran Sasso and operated at ~ 10 mK. A complete description 
can be found in [5] and references therein for what concerns 
the electronics and the DAQ and in l6[] for what concerns the 
cryogenic setup and the shields. 

3. Data analysis 

CCVR data processing, from raw data to the final spectra, 
follows the procedure described in details in 01 . 

The pulse amplitude is estimated by means of an Optimum 
Filter (OF) technique JH] that reduces the noise superimposed 
to the signal, maximizing the signal to noise ratio. 



2 



The mean detector response and the noise power spectral^ 
density, needed to build the OF transfer function, are esti-i; 
mated from bolometric pulses and noise samples (data samplesi; 
recorded randomly and without triggered events) by a properi: 
averaging procedure. « 

Gain instability corrections [7[] are performed using thei: 
5407.5 keV a line from 210 Po. This element is always presents 
in recently grown crystals but decays away with a half life ofic 
138.38 days. n 

The energy calibration is performed using 232 Th y sources 
inserted inside the cryostat external lead shield. An example of l: 
calibration spectrum (Channel2 - CCVR1) is shown in Fig.Q] l: 
Gamma lines from the 232 Th decay chain are clearly visible in 1! 
the spectrum. The calibration is performed using a third-order 11 
polynomial function and the 210 Po peak is used in addition to 1 ' 
the 232 Th y peaks. 




142 
143 
144 
145 
146 
147 



Figure 1: Calibration spectrum of Channel2 - CCVR1. Gamma lines from the 
232 Th decay chain are visible in the spectrum. 

100 

101 Each CCVR sum spectrum is composed by events which sur- 

102 vive two different type of data selection: global and event-based 

103 requirements. 

104 Global requirements are applied following criteria decided a 

105 priori on the detector performances (an excessive noise level, 
loo ADC saturation, etc.). They identify bad time intervals that 
107 need to be discarded. These kind of cuts introduce a dead time 
ioe that is accounted for by properly reducing the live time of the 
loo interested detector. 

no Event-based requirements comprise: pile-up rejection, pulse- 
in shape and coincidence selection. 

112 The presence of a pile-up prevents the OF algorithm from 

113 providing a correct evaluation of the pulse amplitude. The pile- 

114 up rejection is performed by imposing an extendable (paralyz- 

115 able) dead window of 7 seconds to each event. 

no The pulse-shape analysis is used to reject non-physical 
117 events. The pulse shape parameters are the rise time and de- 
ns cay time of the OF-filtered waveform and other parameters that 
no measure the deviation of filtered raw signal from the average 

120 detector response. u 

121 As a first step, each CCVR spectrum is corrected for the cor- 

122 responding efficiency of the event-based cuts (from TableO and 1 ' 

123 then the spectra are summed together. Four types of spectra arei> 

124 produced for all CCVRs: K 

1E 

125 • Total Energy spectrum (TOT): it contains all the general 

126 cuts and the pulse shape cuts. u 



• Anti-coincidence Energy spectrum (Ml): it contains the 
events that caused an energy deposition in one crystal only 
(anti-coincidence cut). For what concerns a particles, this 
corresponds to bulk events (that is decays where the emit- 
ted a particle is absorbed within the crystal itself) and to 
surface events, generated by decays occurring on the sur- 
face of an inert material that hit a facing crystal, or to a de- 
cays occurring on a crystal surface, whose escaping prod- 
ucts are absorbed by inert materials. 

• Coincidence Energy spectrum (M2): it contains the 
events that caused an energy deposition in two crystals. 
For what concerns a particles, these events arise from a 
decays occurring on a crystal surface faced to another crys- 
tal. 

• Coincidence Sum Energy spectrum (M2sum): it con- 
tains the sum energy of multiplicity 2 events. For instance 
in the case of an a decay on the surface of a crystal, both 
the energy of the a particle E\ in the facing crystal and 
the corresponding nuclear recoil energy Ei in the original 
crystal are detected. In the spectrum M2sum the variable 
Ejot — E\ + E2 is plotted. 

Fig. [2] shows the sum spectra of all CCVRs. 




Anti coincidence - Ml 
Anti-Coincidence - M2 



hllibiilhill/w 4 




6000 7000 8000 
TotalEnergy [keV] 



Figure 2: Energy spectra for the full CCVR statistics. Top: anti-coincidence 
and coincidence spectra. Bottom: Sum energy spectrum of multiplicity 2 
events. 



3.1. Efficiency of event-based cuts 

Due to the high rate of 210 Po events, there is a significant 
loss of efficiency due to pile-up rejection. The efficiency is es- 
timated as: 



^pile-up 



l-P 



pile -up 



Jrr-T) 



(1) 



3 



154 where P is the probability of a pile-up, r is the counting rate ofigg 

155 the events that passed the global cuts described above, and T is2oo 

156 the length of a time interval containing an event during whicli2oi 

157 the occurrence of another event would be considered pile-up. 202 

158 The interval T contains some time after the event during which 

159 a double pulse would result and some time before the event 203 

160 during which the event's baseline would be spoiled by the tail 204 

161 of the preceding pulse. Properly, T depends on the energy of 205 

162 the other event: the higher the energy of the other event, the 206 

163 longer its tail remains too large. The counting-rate is channel- 

164 dependent and sometimes also time-dependent. In the specific 207 

165 case of 210 Po events, the counting-rate obviously decreases with 208 

166 time due to polonium decay (half-life: 138.38 days). 

167 For the sake of simplicity, an average pile-up reduction ef-209 

168 ficiency for each CCVR is computed. This will apply to all 

169 channels in the full energy range and is calculated using for- 

170 mula (Q} with T — 7 s and using as r the global counting rate 210 
after the general cuts. 

211 



Run 


^pile-up 


SPSA 


sac 


CCVR1 


0.84 +0.01 


0.96±0.01 


0.987±0.003 


CCVR2 


0. 


88 +0.01 


0.98±0.01 


0.982±0.005 


CCVR3 


0. 


89 ±0.01 


0.97±0.01 


0.990±0.002 


CCVR4 


0. 


88 +0.01 


0.94±0.02 


0.987±0.003 


CCVR5 


0. 


89 ±0.01 


0.98±0.01 


0.990±0.003 



Table 3: Efficiencies of event-based cuts for each CCVR. 



171 

172 The pulse shape cuts efficiency epsa is evaluated on the back- 

173 ground peak at 2614.5 keV due to 208 T1, by a simultaneous fit 

174 on both the spectra of accepted and rejected events as detailed 

175 inH. 

176 The same procedure is applied for the evaluation of the anti- 

177 coincidence cut efficiency £ac- In this case, instead of the 

178 2614.5 keV line (which is usually in coincidence with other y 

179 lines), the photopeak at 1460.8 keV due to 40 K is used. Results 

180 are summarized in Table [3] 

181 3.2. Monte Carlo simulations 

182 In order to extract numerical information on the activity 216 

183 or contamination of a given nuclide in CUORE crystals from2i7 

184 CCVRs data, it is necessary to rely on Monte Carlo simula-218 

185 tions, capable of reproducing the main features of the detector 2 i 9 

186 geometry and response. 220 

187 CCVRs simulations are performed with the GEANT4-based22i 

188 code developed by the CUORICINO ad CUORE collaborator^ 

189 and described in [9] and [10]. The simulation takes into ac-223 

190 count the energy resolution and the threshold of each detector.224 

191 For bulk contaminations, an homogeneous distribution of the225 

192 given nuclide inside the crystals is assumed. For surface con-226 

193 taminations, different depths of the contamination layer (frorri227 

194 0.01 jum to 10 fj.m) are considered and the contamination den-228 

195 sity profile is assumed to decrease exponentially with the depth229 

196 from the surface to the bulk of the crystal. 230 

197 In CCVR1 there is a not negligible number of measurementS23i 

198 (~ 4%) with only 2 active facing crystals. This affects the effi-232 



ciency calculation of surface contaminations. To take into ac- 
count this variation from the standard configuration with 4 crys- 
tals, 2 different efficiencies are calculated for each Monte Carlo 
simulation: 

• e 2 mc'- corresponding to the CCVR1 configuration with only 
2 active crystals; 

• s 4 ^: corresponding to the standard CCVR configuration, 
with 4 active crystals; 

The average efficiency, weighed on the proper lifetime At 1 , is 
computed using the formula: 

pick a flch 1 „4c/j a fAch 

EMC ~ Af 2f/l + Af 4 ^ (2) 

4. Results on background rates of CUORE crystals 

From the energy spectra of all CCVRs the background rates 
in various energy regions can be calculated. 

Six energy regions of interest are identified in the spectra and 
the corresponding count rates for anticoincidence (Ml) and co- 
incidence (M2) spectra are calculated. Results for the global 
spectra are reported in Table [4] (errors are statistical only). 





Continuum 


190p t 


Continuum 




(2700, 3200) 


(3200, 3400) 


(3400, 3900) 




[keV] 


[keV] 


[keV] 


Ml 


0.19+0.02 


0.38+0.04 


0.09+0.01 


Ml-PoSub 


0.13+0.02 


0.34+0.04 


0.06+0.01 


M2 


0.05+0.01 


0.02+0.01 


0.025+0.006 


M2-PoSub 


0.02+0.01 


0.01+0.01 


0.008+0.008 




U/Th 


210p o 


U/Th 




(4000, 5000) 


(5000, 6000) 


(6000, 8000) 




[keV] 


[keV] 


[keV] 


Ml 


0.19+0.01 




0.057+0.004 


Ml-PoSub 


0.13+0.01 




0.057+0.004 


M2 


0.04+0.01 




0.014+0.002 


M2-PoSub 


0.014+0.007 




0.014+0.002 



Table 4: Count rates measured in [counts/keV/kg/y]. Errors are statistical. 



The continuum region (2700, 3200) keV is of great interest 
since it is the region immediately above the Q-value of the neu- 
trinoless double beta decay of 130 Te J3]. 

In the region (3200, 3400) keV a contribution of the a line 
from 190 Pt is expected. This contamination is almost unavoid- 
able for Te02 crystals, as explained in Sec. [5] 

From 4000 to 8000 keV the contribution of the various a 
lines from U and Th decay chains is expected. 

In between, there is the region (5000, 6000) keV, which is 
affected by the 210 Po contamination (see Sec. 15. 11 1. This pro- 
duces not only a peak at the energy of the Po a line 5407.5 
keV but also a broad background over the entire region due to 
mis-identified pile-up events (above the peak energy) or to the 
escape of the a that releases part of its energy in a inert material 
(below the peak energy). An indication of the rate in this region 
is of no particular interest, also because of the relative short 



4 



272 
273 
274 

275 
276 
277 
278 
279 
280 



half-life of 210 Po (138.38 days) that guarantees a huge reduc- 
tion of this count rate when CUORE will start the data taking. 

4.1. Subtraction of 2m Po induced rate 

Because of the presence of 210 Po, an excess of count rate 
could arise in the Ml and M2 spectra below the energy of the 
5407.5 keV a line, if the contamination of 210 Po is close enough 
to the surface for the a to escape and release part of its energy 
in an inert material (Ml spectrum) or in a nearby detector (M2 
spectrum). This contribution from the M2 spectrum can be esti- 
mated, calculating for each energy region the rate of M2 events 
in which the total energy Ejot lies in the interval (5407.5+50) 
keV. The M2 count rate subtracted for this contribution is de- 
fined as M2-PoSub and reported in Table|4]for comparison. 

In a similar way the contribution of surface 210 Po in the Ml 
spectrum is evaluated. Because of the geometry of the CCVR 
setup, our coincidence analysis is sensitive only to 8 faces over 
24. So the anti-coincidence count rate due to 210 Po of the re- 
maining 16 faces is assumed to be two times the coincidence 
one, calculated on the M2 spectrum as described above@ The 
Ml count rate subtracted for the surface 210 Po contribution is 
defined as Ml-PoSub and reported in Table [4] 

4.2. Background rates comparison with previous TeC>2 detec- 
tors 282 

It is interesting to compare CCVRs rates with the ones ob-283 
tained with previous detectors, like TTT ifTHl or Cuoricino J3l. 
This comparison is shown in Table [5] where the Ml and M2 
count rates for CCVRs are after polonium subtraction, as de 
scribed in 14.11 

285 

The notation (2700, 3900) refers to the combination of both 
the continuum region (2700, 3200) keV and (3400, 3900) keV, 286 
excluding the 190 Pt energy region of (3200, 3400) keV. It can be 28 7 
inferred that : 

288 

• in the region (2700, 3900) keV the CCVRs anti- 
coincidence rate is compatible within 1.8 cr with the cor- 
responding TTT value. 290 

• in the region (4000, 5000) keV, as already measured in 291 
the TTT run, a reduction in the count rate with respect to 292 
Cuoricino is observed. 293 

• in the region (5000, 6000) keV a comparison is not possi- 294 
ble since CCVR is affected by the high rate of 210 Po. This 295 
is due to the fact that in CCVRs only recently grown crys- 296 
tals are measured, unlike TTT or Cuoricino. J 

298 

• in the region (6000, 8000) keV, as already measured in 299 
the TTT run, a reduction in the count rate with respect to 30 o 
Cuoricino is present. However, the CCVR count rate is 301 
greater than the TTT value, probably because of the pres- 302 
ence of mis-identified pile-up (in Ml) or coincidences (in 
M2) with 210 Po events, extending above 6000 keV. This 303 
contribution should decay away with 210 Po. 304 
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210n 


U/Th 






(5000, 6000) 


(6000, 8000) 






|keV J 


71 

[kev J 


CCVR 


Ml 




0.057+0.004 




M2 




0.014+0.002 


TTT 


Ml 


1.30+0.07 


0.025 + 0.004 




M2 


0.09±0.01 


0.005 + 0.002 


Cuoricino 


Ml 


0.846+0.004 


0.099 + 0.001 




M2 


0.173+0.002 


0.0163 ± 0.0004 



Table 5: Count rate comparison with previous detectors, measured in 
[counts/keV/kg/y]. Here CCVR values are after polonium subtraction, as de- 
scribed in 14.11 Errors are statistical. 



A more detailed evaluation of the comparison can be found 
ind. 



In some of CCVR1 measurements (~ 4% of the CCVR lifetime) only two 
facing crystals were active. The extrapolation of M2 counts to Ml is performed 307 
considering that the coincidence analysis is sensitive only to 2 over 12 faces. 308 



284 5. Results on bulk contaminations of CUORE crystals 



Bulk contaminations in CUORE crystals are expected from: 

• 210 Po, as a result of the chemical affinity between polo- 
nium and tellurium; 

• natural contaminants, like 238 U and 232 Th with their ra- 
dioactive chains; 

• 210 Pb, as a result of the deposition of lead nuclei produced 
by 222 Rn decays, during the crystal handling in free atmo- 
sphere; 

• 190 Pt, due to the fact that platinum is used in several phases 
of the crystal production cycle. Platinum crucibles are 
used for the calcination of Te02 powder used for the crys- 
tal growth and the growth crucibles are made of platinum 
foil. The central part of the as-grown crystal ingot is se- 
lected for the CUORE crystals in order to avoid the risk of 
platinum contamination on the surface of the crystal due 
to possible diffusion during the growth process. Possible 
Pt contaminations are therefore in the bulk of the CUORE 
crystals. 

For 210 Po the activity is determined from a fit to the anti- 
coincidence rate vs. time behaviour (see Sec. 15. U . 

For 238 U and 232 Th (see Sec. I5.2t . the limit on the level of 
contamination is determined from the intensities of the a peaks 
in the anti-coincidence spectrum of all crystals, or from inte- 
grals centered at the peak position if the peaks are not visible. 



5 



This is because an a decay from bulk contamination releases 
the entire Q-value of the reaction (a energy + nuclear recoil) in 
a single crystal. 

For 210 Pb, the limit on bulk contamination is determined 
from a fit in the energy region (40, 60) keV in a subset of 
CCVRs data with high statistics and low threshold (see Sec. |6). 

5.1. 210 Po bulk activity 

The activity of 210 Po can be measured from the intensity of 
5407.5 keV a line in the anti-coincidence spectrum. 

The plot in Fig. [3]shows the global rate of 210 Po events over 
time for CCVR1 . For each channel the 210 Po events are selected 
in a ± 20 keV window around the energy of the a line. 
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^ 6000 
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Figure 3: Global rate of 210 Po events over time for each CCVR1. The fit result 
is overlayed. 

The units on the x-axis are days since the start of the first 
background measurement. Each point represents a group of 
measurements whose livetime is at least 5 days. 

The horizontal error bars indicate the beginning and the end 
of each group of runs, they are for visualization only. 

Each point is corrected with the corresponding rate-based ef- 
ficiency, calculated with the equation (Q]i. A larger dead time 
window (9 seconds instead of 7 seconds) compared to the rest 
of the analysis is used, to be more conservative in the removal 
of pile-up pulses. 

The fit function is a pure exponential: 

J-ln2 t/T m ) 



r(t) = r e [ 



(3) 



where ro is the rate at the beginning of the measurement and 
T\/2 is the 210 Po half-life. The half-life of the exponential decay 
has been evaluated for all CCVRs and it is shown in Tab[6j it is 
in good agreement (1 tr) with the half-life of 210 Po (138.38 d). 3! 
This indicates that the 21() Po contamination is out of equilibrium 
and it is not being fed by 21l) Pb. 31 
From the value of ro returned by the fit for each CCVR and 
for each crystal, the 210 Po activity at the beginning of the mea- 
surement is extracted as follows: 

3! 

ro 



A [Bq/kg] = 



(4> 



86400[s/day] m[kg] 

where m is the crystal mass. 3e 
Knowing the time elapsed since the "crystal birth date"3E 
(growth completed, before the cut and shape) and the start ofiiE 
the measurement, the 210 Po activity at production is computed. 3e 



Pun 


Half lifp 




[days] 


CCVRl 


138.0±2.5 


CCVR2 


132.3 + 14.7 


CCVR3 


139.4 ±6.3 


CCVR4 


137.8 + 11.5 


CCVR5 


136.6 ±9.1 



Table 6: 210 Po half-life for all CCVRs. The activity has been fitted using an 
exponential function. All the values are consistent (within 1 <x) with 138.38 
days. 



The 210 Po activities for each CCVR and for all crystals are 
reported in Table Q The results for crystals 007 and 011, mea- 
sured first in CCVRl and again in CCVR2, are consistent. 

All crystals tested in the first five CCVRs are well below the 
limit of 0.1 Bq/kg imposed to the crystal producers. 



CCVR 


Crystal 


21u Po activity [Bq/kg] 


1 


041 


0.0257 ± 0.0001 


1 


011 


0.0510 ±0.0005 


1 


039 


0.0229 ± 0.0001 


1 


007 


0.0414 ± 0.0004 


2 


076 


0.021 ± 0.004 


2 


011 


0.07 ± 0.02 


2 


096 


0.055 ± 0.006 


2 


007 


0.047 ± 0.015 


3 


190 


0.0078 ± 0.0005 


3 


236 


0.0159 ± 0.0005 


3 


180 


0.0203 ± 0.0008 


3 


229 


0.0283 ± 0.0008 


4 


340 


0.032 ± 0.004 


4 


313 


0.005 ± 0.001 


4 


354 


0.039 ± 0.004 


4 


380 


0.040 ± 0.004 


5 


455 


0.019 ± 0.002 


5 


416 


0.024 ± 0.003 


5 


436 


0.032 ± 0.003 


5 


421 


0.020 ± 0.003 



Table 7: 210 Po activity at production time for all CCVRs crystals. 



5.2. U/Th bulk contaminations 

238 U and 232 Th bulk contaminations are evaluated from the 
anti-coincidence sum spectrum of all CCVRs. For each of the 
peaks reported in Table [8] the number of counts (corrected by 
the efficiency of the event-based cuts from Table [3} within an 
energy window of ±6 <x around the a Q-value is estimated. A 
cr of 2.2 keV, corresponding to the average between the values 
of the 210 Po peak <xs for each CCVR, weighted by the corre- 
sponding lifetime of that run (see Table [9] for details) is used. 
Using the Bayesian approach, the upper limits N„ at 90% C.L. 



6 



v^iidin 


Nuclide 


Energy 
[keV] 


odii-iiie 


AT 

90% C.L. 




l_.fl din 


Nuclide 


Upper limit 
[Bq/kg] 


Upper limit 

[g/g] 


238t t 


238tt 


4270.0 


4.47E+09 y 


11.7 






238i t 


2.5E-07 


2.0E-14 




234 u 


4858.8 


2.45E+05 y 


20.8 






234 u 


4.7E-07 


3.6E-14 




230 Th 


4770.0 


7.54E+04y 


24.9 






230 Th 


5.7E-07 


4.4E-14 




226 Ra 


4870.6 


1599 y 


30.0 






226 Ra 


6.7E-07 


5.3E-14 




218p 


6114.7 


3.05 min 


7.4 






218p 


1.6E-07 


1.3E-14 


™Th 


232 Th 
212 Bi 


4082.8 
6207.1 


1.4E+10y 
60.55 min 


5.9 
13.5 




™Th 


232 Th 
212 Bi 


1.3E-07 
8.4E-07 


3.1E-14 
2.1E-13 



Table 8: 90% C.L. limits on the number of events ascribed to several nuclides 
from uranium and thorium decay chain. For each nuclide is also shown the 
Q-value and the half-life of the a decay. 



Table 10: Upper limits at 90% C.L. on the activity and on the bulk contamina- 
tion of uranium and thorium decay chains in the hypothesis of secular equilib- 
rium. 



Run 


Livetime 


<T 




[days] 


[keV] 


CCVR1 


59.9 


1.6 


CCVR2 


19.4 


1.4 


CCVR3 


43.05 


2.2 


CCVR4 


25.8 


3.5 


CCVR5 


30.3 


2.8 



386 
387 
388 
389 
390 



Table 9: Summary of CCVRs energy resolutions, evaluated with a gaussian fit 
of the 5407.5 keV peak from 210 Po. The average value of 2.2 keV, used for 392 
the global analysis is obtained averaging the resolution of each CCVR with the 
corresponding livetime. 



362 are obtained, assuming expected background counts and a flats^ 

363 prior for the signal 11211 (see TablelSl. 397 

364 The upper limit on the activity for each nuclide is calculated 393 

365 using the following formula: 399 



A u [Bq/kg] = 



e M c T[s] m[kg] T 



(5)* 



where smc is the Monte Carlo detection efficiency, T the 403 
CCVRs livetime, m the crystal mass and T the branching ratio 404 
of the nuclide. In this analysis a total containment of the anti- 405 
coincidence events in the crystals is assumed (smc =!)■ Results 406 
are shown in Table [TOl 407 

The upper limit for confidence level of 90% for U/Th bulk 408 
contaminations are then calculated in the hypothesis of secular 409 
equilibrium within the uranium and thorium decay chains. Re- 410 
suits for the total CCVRs data are shown in the last column of 411 
Table [TO] 

In the most conservative approach, the bulk contamination 413 
limit on 238 U and 232 Th is set considering the most active nu- 414 
elide for each chain. The upper limit at 90% C.L. for uranium 415 
and thorium bulk contamination are: 416 



23N 



u 



232 



Th 



< 5.3 

< 2.1 



10- 14 [g/g] 
10- 13 [g/g] 



Both values of the upper limits are within the contract specifica-421 
tion of 3T0~ 13 [g/g]. For the 238 U decay chain, the contribution422 
from 210 Pb is treated separately (see Sec.|6]i. 423 
It is important to stress the fact that the above upper limitS424 
are calculated under the hypothesis that the observed counts for425 



each nuclide are entirely due to a bulk contamination of that 
nuclide. This is a conservative hypothesis since there is not a 
clear indication (for example a line) that such a bulk contami- 
nation actually exists and that the observed counts are not due 
to background of some other origin. 

6. 210 Pb activity 

During the production of CUORE crystals great care is de- 
voted in order to minimize the exposure of the crystals to free 
atmosphere, to avoid contamination of radon and its daugh- 
ters. Nevertheless, a small 21() Pb contamination (bulk or sur- 
face) might occur. This contamination can be very dangerous 
because, even if the soft beta (16.96 keV and 63.5 keV) and 
y (46.5 keV) radiation from 210 Pb is generally self absorbed, 
the induced bremsstrahlung of the daughter 210 Bi can give rise 
to a continuum up to 1.16 MeV Moreover, the a decay of the 
daughter 210 Po can contribute to the continuum background in 
the double beta decay energy region. 

In CCVRs crystals the 21() Pb contamination cannot be esti- 
mated from 2I0 Po, because this contamination is out of equi- 
librium (see Sec. 15. U . The only available signature is a com- 
bination of a beta spectrum (end point 16.96 keV) with a de- 
excitation energy of 46.5 keV (most of the cases through a con- 
version electron). This produces a broad signature in the energy 
region (40, 60) keV, whose shape depends on the location of the 
Pb contamination. 

Several simulations are performed both for a bulk contami- 
nation and for a surface contamination with exponential density 
profile [9] and contamination depth varying from 0.01 /zm to 
10 fim. Even if, for the sake of simplicity, a common energy 
threshold of 50 keV is set for all CCVRs, there are subset of 
data where the threshold can be set to a lower value in order to 
look for this signature. CCVR1 data are used because they struk 
a balance between high livetime (CCVR1 has the highest statis- 
tic) and good bolometer performances. Only for CCVR1 data, 
the analysis is repeated with an energy threshold of 40 keV. 
The anti-coincidence spectrum between 40 and 60 keV is fit- 
ted with an exponentially decreasing background and allowing 
the presence of a 21() Pb spectrum with shape taken from Monte 
Carlo simulation. The free parameters of the fit are the two pa- 
rameters of the exponential background and the total number 
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Figure 4: Top: Monte Carlo simulation of a 210 Pb bulk contamination in 
CCVR1. Bottom: CCVR1 low energy anti-coincidence spectrum. The energy 
threshold is set to 40 keV. The red line represents the fit to the spectrum with 
an exponential background and allowing the presence of a 210 Pb bulk contami- 
nation with the above shape. No hint of such contamination is found. 



426 of counts from 210 Pb. The fit is repeated for each Monte Carlo 

427 signature. The number of 210 Pb events is compatible with zero 

428 within the error for all signatures. An example of Monte Carlo 

429 spectrum for a 210 Pb bulk contamination and the corresponding 

430 fit to the experimental spectrum is shown in Fig. [4] 

431 The upper limit at 90% C.L. on the number of 210 Pb counts 

432 is N u = 1 .644 cr, where <x is the error on the number of counts 

433 returned from the fit. The upper limit on bulk contamination is 

434 computed using Eq. [5] The upper limit on surface contamina- 460 

435 tion is extracted using the following formula: 

462 

A„[Bq/cm 2 ] = Ntt (6) 463 

s M c T[s] 5 [cm 2 ] 464 

437 where smc is the Monte Carlo detection efficiency, T is the live- 4 ^ 

438 time and S is the surface of a crystal (150 cm 2 ). Results are 466 

439 shown in Table QT| The upper limit on the activity for bulk«67 

440 contamination is below the contract limit of 10~ 5 Bq/kg (see*™ 

441 Table Q]). 469 

470 

442 7. Results on surface contaminations of CUORE crystals 471 

472 

443 Surface contaminations of CUORE crystals are expected473 

444 from the same nuclides listed in Sec. [5] except for 190 Pt, which474 

445 is all in the bulk. Uranium and thorium surface contaminations475 

446 can be investigated by means of coincident events in two fac-476 

447 ing crystals: surface contamination of an ff-decaying nuclide477 

448 should appear in the total energy spectrum (M2sum in Fig. |2])478 

449 as a peak at the Q-value of the decay because the total energy (am 

450 + nuclear recoil) is collected by the two facing crystals. At the4so 

451 same time, if the surface contamination is deep enough, the a4si 

452 particle can be absorbed by the crystal, giving rise to the same482 

453 signature of a bulk event. 483 



contamination 


Upper limit yu% 


bulk 


3.3E-06 


surf. 0.01//m 


9.8E-07 


surf. 0.1//m 


3.8E-08 


surf. 0.2//m 


2.2E-08 


surf, lyt/m 


9.2E-09 


surf. 5pm 


5.6E-09 


surf lOyum 


4.9E-09 



Table 1 1 : Upper limits on the activity of bulk and surface contamination of 
210 Pb. Bulk contamination is given in [Bq/kg], surface contaminations are 
given in [Bq/cm 2 ]. 



454 7.1. U/Th surface contaminations 



The M2sum spectrum in Fig.[2]shows a huge background for 
energies above the polonium a line due to random coincidences 
between 2I() Po and low energy events. At a closer inspection 
(see Fig. [5j, the peaks originated by the sum of a 5407.5 keV a 
line and a low energy y line arising from Te metastable isotopes 
are clearly visible. See also Table [T2lfor details. 



Te isotope 


Energy 


Energy + 5407.5 




[keV] 


[keV] 


m Te 


88.3 


5495.8 


129 Te 


105.5 


5513 


125 Te 


144.8 


5552.3 


123 Te 


247.5 


5655 


121 Te 


294 


5701.5 



Table 12: Energy of the gamma lines from Tellurium metastable isotopes and 
corresponding sum energy when in coincidence with a 210 Po event. 



Because of the presence of these coincidences, only nuclides 
with a Q-value lower than 5407.5 keV (see Table [TBI are used 
for surface contamination analysis. 2I8 Po and 212 Bi are there- 
fore discared. The scatter plot of M2 events is shown in Fig. [6] 
The shadowed region contains the events with a total energy 
within 4 and 5 MeV. 

For each nuclide listed in Table [T3l the number of counts is 
computed from the M2sum energy spectrum in an energy win- 
dow of + 6cr around the Q-value and divided by s 2 (where s 
is the efficiency of the event-based cuts of Table [3}, since two 
coincident events have independent probabilities of passing the 
cuts. In Table [T3l are shown the corresponding upper limits at 
90% C.L., computed using the Bayesian approach with ex- 
pected background counts and a flat prior for the signal lfl2ll . 

Monte Carlo spectra for 238 U and 232 Th contamination on 
the crystal surface are generated with exponential profile and 
various penetration lengths. For each contamination depth, the 
containment efficiency is calculated in a + 6<x interval around 
the Q-value of each nuclide both for the M2sum spectrum and 
for the Ml spectrum. Results are shown in Fig. [7] For penetra- 
tion length of 0.01 and 0. 1/j.m the higher containment efficiency 
comes from the M2sum spectrum, whereas for depths of 1, 5 
and 10 pm this arises from the Ml spectrum. The efficiency 
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Figure 5: Top: M2sum spectrum in the region (5300, 6000) keV. Peaks origi- 489 
nated by the coincidence of the 5407.5 keV a line and a low energy y line from4 90 
Te metastable isotopes are clearly visible. Bottom: total and anti-coincidence 
spectra in the region 0-500 keV: low energy lines from Te metastable isotopes 
are clearly visible. 492 




4000 5000 6000 

TotalEnergy-Energy [keV] 



Figure 6: Scatter plot of M2 events. The shadowed region contains events 
total energy between 4 and 5 MeV, used for surface analysis. 
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495 
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498 
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507 
508 



with 



Chain 


Nuclide 


Half-life 


Energy range 


N u 






[y] 


[keV] 


90% C.L. 


m v 




4.47E+09 


4257 - 4283 


4.6 




234 u 


2.45E+05 


4845-4871 


4.6 




230 Th 


7.54E+04 


4757 - 4783 


4.6 




226 Ra 


1599 


4857 - 4883 


10.5 






1.4E+10 


4069 - 4095 


2.3 



Table 13: Upper limits at 90% C.L. on the number of counts ascribed to several 517 
nuclides from uranium and thorium decay chain from the M2sum spectrum. 518 
For each nuclide, the half-life of the a decay is also shown. 



Th232 Ml 
Th232 M2sum 
U238 M1 
U238 M2sum 



10 

Depth [micronj 



Figure 7: For each penetration length, the containment efficiency of surface 
events is computed for both Ml and M2sum Monte Carlo spectra. 



containment value for 0.2 fim depth is very similar both for Ml 
and M2sum spectrum. 

The drop in Ml containment efficiency for very thin layers 
of contaminations is due to the fact that, of the about 50% of 
decays where the alpha is emitted toward the bulk of the crystals 
and the recoil is emitted towards outside, in most of the cases 
the recoil escapes undetected, thus the full Q-value of the alpha 
decay is not recorded. 

The upper limits at 90% C.L. for the surface activity of each 
nuclide are evaluated using Eq.|6]where N u is the 90% C.L. up- 
per limit on the number of observed events from the M2sum 
or the Ml spectrum (see Table [131 for the M2sum counts and 
Table [8] for the Ml counts) and smc is the Monte Carlo aver- 
age efficiency defined in eq. [2] for the corresponding spectrum 
(M2sum or Ml) and for the given signature. The confidence 
intervals for surface activity contaminations are calculated for 
both the M2sum and Ml spectra, normalizing each signature 
with the corresponding Monte Carlo efficiency. The signature 
giving the most stringent result is taken into account. 

For surface contaminations of 0.01, 0.1 and 0.2 fim depths, 
the most stringent limits come from the M2sum spectrum. For 
the remaining depths (1,5 and 10 pun) the surface activity re- 
duces practically to a bulk activity and the Ml signature pro- 
duces the most stringent limits (see Table[8]for details). 

For the 238 U chain, the surface contamination for 2 peaks is 
evaluated: 

• 238 U, the chain parent; 

• 226 Ra, the most active line both in Ml and M2sum spectra. 

The contribution from 210 Pb is treated separately (see Sec. [6j. 

For the 232 Th chain there is only one useable line for this 
analysis, that is the one from 232 Th. This means that there is 
no way of testing the portion of the chain below 220 Rn and take 
into account a possible non-equilibrium of the chain, as done 
for the bulk contamination. 

The results for surface contaminations are shown in Tablefl4l 
As explained in Sec. 15. 21 the above confidence intervals are cal- 
culated under the hypothesis that the observed counts for each 
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544 
545 



nuclide are entirely due to a surface contamination of that nu- 
clide in the corresponding depth. Again, this is a very conser- 
vative hypothesis. 



Depth 


Nuclide 


Upper limit 90% C.L. 






[Bq/cm 2 ] 


O.Olyum 


a»U 


3.1E-09 




226 Ra 


6.3E-09 




232 Th 


1.6E-09 


0.2yum 


mjj 


3.8E-09 




226 Ra 


7.6E-09 




232 Th 


2.0E-09 




238tj 


3.7E-09 




226 Ra 


8.9E-09 




232 Th 


1.9E-09 


5/jm 


238JJ 


2.0E-09 




226 Ra 


5.4E-09 




232 Th 


1.0E-09 


lOyum 


238JJ 


1.7E-09 




226 Ra 


4.4E-09 




232 Th 


8.3E-10 



552 
553 
554 
555 



Table 14: Upper limits at 90% C.L. for surface contamination, for different 
penetration length values. See text for details on the calculation of confidence 
intervals. 557 



559 

560 



524 8. Extrapolation to CUORE background 

525 In order to evaluate the contribution to the CUORE back- 

526 ground arising from crystal contamination, a simulation both 

527 for bulk and surface contamination reproducing the CUORE 

528 geometry and studying the contribution to the double beta de- 

529 cay (DBD) energy region is performed. 

530 Being interested in a conservative upper limit to the CUORE 

531 background, in this extrapolation it has been assumed that 

532 CUORE crystals will have the same activity as crystals tested 

533 in CCVRs. However, when CUORE will start the data taking, 

534 most of the 21() Po (and other short lived states) will have de- 

535 cayed and the efficiency for pile-up rejection will be increased. 

536 8.1. Background from bulk contamination 

537 The CUORE geometrical efficiency for a uniform bulk con- 

538 tamination of crystals in 21() Pb, 238 U and 232 Th is estimated via 

539 Monte Carlo, and the corresponding CUORE background at the 

540 DBD energy region (Q-value + 30 keV) is calculated using the 

541 formula: 

a <XUORE-bulk 
AcCVR-bulk S MC 



bkg 



CUORE 



AE 



(7) 



where AccvR-tmik are the values of 210 Pb, 238 U and 232 Th ac- 
tivities from Tables [10] and QT] e cuoRE-buik is eva i uate( j f rom 
Monte Carlo simulation and AE = 60 keV. 

The results are shown in Table [151 
In the most conservative approach, considering the most ac-sei 
tive line ( 212 Bi), the following upper limit to the CUORE back-562 
ground at the DBD energy due to bulk contamination of crystalss63 
is set: 1.1 • 10~ 4 counts/keV/kg/y. 564 



chain 


Nuclide 


Upper Limit 90% C.L. 
[counts/keV/kg/y] 




2io pb 


2.6E-05 


238TJ 


238TJ 


8.1E-07 




234TJ 


1.4E-06 




230 Th 


1.7E-06 




226 Ra 


2.1E-06 




218p 


5.1E-07 


232 Xh 


232 Xh 


1.7E-05 




212 Bi 


1.1E-04 



Table 15: Extrapolation to CUORE background from CCVRs bulk contamina- 
tion limits from Tabs. llOl and ll II 



8.2. Background from surface contamination 

In a similar way, the CUORE geometrical efficiency for a 
surface contamination for several depths is estimated. The cor- 
responding CUORE background is extrapolated using the for- 
mula: 



bkg ( 



iCCVR-surf E 



.CUORE-surf 



CUORE 



AE Mcuore 



(8) 



blesQT]and[T4] e™ 0RE ~ surf is estimated via Monte Carlo, S is 
the surface of the CUORE crystals, AE = 60 keV and Mcuore 
= 0.75 • 988 kg. 



Depth 


Nuclide 


Upper limit 90% C.L. 






[counts/keV/kg/y] 


O.Ol^m 


2i0p b 


1.6E-03 




238TJ 


8.9E-05 




226 Ra 


1.8E-04 




232 Th 


1.8E-05 


O^m 


210p b 


2.3E-04 




238TJ 


3.9E-04 




226 Ra 


7.7E-04 




232 Xh 


1.5E-04 


lyum 


2l(.p b 


3.8E-04 




238TJ 


1.2E-03 




226 Ra 


2.9E-03 




232 Th 


4.4E-04 


5yum 


2l0p b 


4.3E-04 




238TJ 


1.3E-03 




226 Ra 


3.4E-03 




232 Th 


4.7E-04 


10/zm 


210p b 


2.9E-04 




238TJ 


7.1E-04 




226 Ra 


1.9E-03 




232 Xh 


3.0E-04 



Table 16: Extrapolation to CUORE background from CCVRs surface contam- 
ination limits from Tables fTTI and 1141 



The results are shown in Table [16] Based on these results, 
in the most conservative approach, a crystal surface contam- 
ination in CUORE within the limits estimated from CCVRs 
data, would yield a rate in the DBD energy region lower than 
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Figure 8: For each penetration depth, the upper limits at 90% C.L. for the c 
rate at the DBD energy region are plotted. The values are computed on 2 
line, and are due only to surface crystals contaminations. 
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565 5.5 • 10 counts/keV/kg/y. This value corresponds to the worst 

see case: a contamination in 238 U and 232 Th at 5/im and a contam- 

567 ination in 210 Pb at 0.01/im. This result is strictly dependent on 

sea the density profile assumed for the contaminants distribution, 

569 as can be seen in Fig. [8] where the upper limit at the 90% C.L. 

570 for the rate in the DBD energy region (computed on 226 Ra line) 

571 is plotted in function of contamination depth. 



1.1-10 4 counts/keV/kg/y and 5.5- 10 3 counts/keV/kg/y re- 
spectively. 
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572 Conclusions 

573 The CUORE Crystal Validation Runs (CCVRs) have been 

574 carried out since the end of 2008 at the Gran Sasso National 

575 Laboratories, in order to test the performances and the ra- 

576 diopurity of the TeC<2 crystals produced at SICCAS (Shanghai 

577 Institute of Ceramics, Chinese Academy of Sciences) for the 

578 CUORE experiment. 

579 In this work the results of the first 5 validation runs, are pre- 

580 sented. The items under analysis are the background rate, the 

581 bulk and surface contamination from several nuclides and the 

582 extrapolation to the CUORE background. 

583 The CCVRs background rate shows a reduction with respect 

584 to the one measured in CUORICINO in all the energy regions 
sas considered. The bulk activity of 210 Po is measured to be within 
sae the limit specified in the contract with the crystals producer. 
587 No indication of a bulk contamination from uranium and tho- 
saa rium decay chains, as well as from 210 Pb (out of equilibrium), is 
sag found. The upper limits are calculated to be within the contract 

590 specification. No indication of a surface contamination from 

591 uranium and thorium decay chains, as well as from 21() Pb (out 

592 of equilibrium), is found. Upper limits at 90% C.L. are calcu- 

593 lated for surface contamination from several nuclides and for 

594 different contamination depths. 

595 An extrapolation to CUORE background from both bulk and 

596 surface contaminations is performed, in the most conservative 

597 assumption that the CCVRs observed background is entirely 

598 due to the bulk and the surface contamination respectively. 

599 Considering only the contribution from bulk and surface crys- 
eoo tals contaminations, the following upper limits to the CUORE 

601 background index in the energy region around the Q-value 

602 of the neutrinoless double beta decay of 130 Te are calculated: 
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